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Abstract 

The derivation of the Eliashberg type equations for a superconductor with 
strong electron correlations and electron-phonon interaction has been presented. 
The proper account of short range Coulomb interactions results in strongly 
anisotropic equations. Possible symmetries of the order parameter include s,p 
and d-wave. We found the carrier concentration dependence of the coupling 
constants corresponding to these symmetries. At low hole doping the d-wave 
component is the largest one. 
PACS: 71.27.+a, 74.20.-z, 74.25.Kc, 75.10.Lp 
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I. INTRODUCTION 



Recently there appeared in the literature a number of papers dealing with the problem 
of electron-phonon interaction in strongly correlated superconductors0~0. It turns out that 
interplay between Coulomb and electron-phonon interactions is very subtle and may result in 
an unexpected behaviour of various quantities0~El. 

The main motivation for all these studies come from experiments on high T c superconduc- 
tors which in one or other way have pointed out the importance of electron-phonon interaction. 
There have been found Raman and infrared-active mo des0 with strongly temperature depen- 
dent linewidth. Small, but nonzero, isotope shift exponent a (T c oc M a , where M the ion 
isotope mass)lii@ is the strong indication of the electron-phonon coupling. Similar conclu- 
sions can be inferred from tunneling and photoemissionlll, neutronEl'Ill, specific heatHl, thermal 
conductivity^!! and other experiments. 

It is the purpose of this paper to derive Eliashberg typeil equations valid for the strongly 
correlated superconductor with electron-phonon interaction. To this end we assume the validity 
of Migdafl theorem which means that we assume relatively weak electron-phonon interaction. 
The strong correlations described by the Hubbard model will be treated via slave boson method. 
The derived equations are strongly anisotropic and lead to the possibility of various types of 
symmetries of order parameter. The relative stability of various symmetries does depend on 
the carrier concentration and other parameters. 

In Section II we present the model Hamiltonian for strongly correlated systems and discuss 
the application of slave boson approximation to it in the U = oo limit, both in mean field and 
beyond it. The derivation of the Eliashberg equations on the real frequency axis is described 
in considerable detail in Section III. Section IV contains presentation of results and discussion. 



II. THE MODEL 

We start with the one band Hubbard model with general form of electron-phonon interaction 
term 

ija ijscr i 
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Here cf a (ci a ) is the creation (annihilation) operator for a spin a electron at site % of the lattice, 
is the a-th component of the displacement vector of the ion, h ia = cf^c^, \i denotes chemical 
potential, is the hopping integral assumed to take on nonzero value —t for i,j being nearest 
neighbour sites, U is the Hubbard on-site repulsion of carriers. if p h denotes the Hamiltonian 
of the lattice. The electron-lattice interaction described by the second term in (1) has two 
components 

1 ijs ~ 1 ijs i v js°iji 

the first of which, T°j s has been derived from modulation of the hopping integra tij in the 

deformed lattice. It is related to the derivative of taken at equilibrium position of an ion 

dt- dt- 
ijs dR% is dRf° js - 

The second part of the interaction i.e. the term V^ s is connected with fluctuations of the crystal 

field. Due to ionic character of high temperature superconductors this term is expected to be 

more important and has to be taken into account even if we have assumed the equilibrium 

value of crystal field = ta to be zero. The systematic derivation and the discussion of the 

electron-ion interaction in context of superconducting oxides can be found in Ref.il. 

We are interested in the strong correlation limit, characterised by the large U values (U ^> 

W, where W = 8t is the width of electron band in a two dimensional square lattice). In 

this limit and for less than half filled band (i.e. concentration of carriers n < 1) it is very 

inconvenient for two carriers to occupy the same site i. It means that double occupation of 

sites is prohibited and the dynamics of carriers is limited to empty and singly occupied states. 

Thus we have to project out all the doubly occupied states of the system. This procedure is 

most conveniently carried out with the help of auxiliary particles§l and in the U = oo limit. One 

rewrites the electron operators in (1) in terms of new fermion operators cf a (cj CT ) and auxiliary 

boson operator bf(bi) as cf a dj a — > cf^Cjabibj' . The term U may then be dropped out at the 

expense of introducing at each site a constraint 

£ ctc la + bfh = 1 (2) 

via Lagrange multiplier A,. The constraint allows for at most single occupation of each size. 
The Hamiltonian (|l|) then becomes 

H = H UjC+c^bibj - fiJ2 c td* + £ JfX^W/ 

ijcr id ijsaoL 

+ E + E A *(E 4*a« + Hh - 1) + #ph- (3) 



Note that on-site terms are not modified. It is connected with the meaning of the boson 
operators bf as operators creating empty states at site i. The hopping of a real electron from 
occupied state j to an empty state at i, described by term cf a dj a in ([!]) consists of motion of 
a fermion between these sites (cf a Cj a ) and at the same time the empty site moves from site i to 

3 

It has been shown that mean field description can be obtained by assuming average values 
(bi) = (bj~) = r and Aj = A at each site. Two parameters r and A entering the mean field 
Hamiltonian can be choosen so as to make minimal the ground state energy Eqs — (H). One 
gets0 

„2 
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"A = i E^fe) + ^ E T^iKcfcj*)- 

ijtj ijsaa 

Here n = 4 Y.ia( c ta c (r) denotes the concentration of electrons in the band (n < 1). There are 
two modifications of the spectrum of electrons encountered on the mean field level. First is the 
band narrowing described by r 2 an its shift described by A. The spectrum of noninteracting 
fermions in the mean field is given by (r 2 tk — /i + A) instead of (e^ — /i) of original electrons (at 
U — 0). For the half filled band n — 1 the system is localised (r = 0). Here is the Fourier 
transform of tij. 

To go beyond mean field approximation for slave bosons it is convenient to define boson 
fluctuating fields via: bf = r + 5bf, bi = r + 5bi. 

In view of our main goal in next section which is derivation of the Eliashberg equations, we 
shall write down the Hamiltonian in terms of Gorkov-Nambu^ field operators: 

V>i = l \ and if>f = (c+pCj-J. 

Using this in (1) we obtain 

H = H MF + H' + H" , (4) 

where in site representation 

ij ijsa 

+ AE(^ + ^ + r 2 -l) + # p \ 

i 

H' = 5>^+t 3 ^ + ASyWbi + 6b+) + *£rT£ s uyt^j(8bi + 5b+) , (5) 

ij ijsa 

H" = AJ2 SbfSk + E h.r; r.r^dbj + E T'^u^fh^ 61*5% . 
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First term of H' describes fermion-boson interaction, while the last one fermion-phonon-boson 
interaction. More complicated interactions are contained in H" . They will be treated in 
a mean field type of approximation. The important point is that H" significantly contributes 
to dynamics of fluctuating boson field. In the mean field approximation for electrons the boson 
fluctuations are described by effective Hamiltonian 

H B = YXASy + tyic+Cj^Sb+Sbi (6) 

i 

For the sake of completeness we write down the Hamiltonian for phonons which we assume to 
be described in the harmonic approximation. In the second quantised form it reads 

H ph = hu qv (a+ v a qu + -) . (7) 

qv * 



III. EQUATIONS OF SUPERCONDUCTIVITY 



To properly describe the superconducting state in the system at hand one has to work in site 
representation. The important point is that in considered U = oo limit the double occupation 
of a given site is strictly forbidden. This means inter alia that correlation functions (c^Cji) 
describing superconducting pairs vanish exactly for i = j, i.e. the on-site pairing is forbidden. 
On the other hand the correlations of the type (c^q ct ) measure the average number of carriers 
at site i, and are allowed to enter into formula. This important fact has first been noted by 
Zieli/nski and coworkers^ and leads, as we shall see, to severe changes in the form of Eliashberg 
equations. 

To derive them we use the equation of motion method for the double time thermodynamic 
Green's functions. For arbitrary operators A an B it reads (we employ here the Zubariev 
notationii) 

w<A|B> w =([A,B]±)+<z:[A,HUB^ u} 
= ([A,B]±)-<A\[B,H\_^> U , 

where <C A\B 3> w is the Fourier transformed, frequency dependent, retarded Green's function. 
[A, B]± denotes anticommutator (+) or commutator (— ). 

Standard procedure^ leads to the following equation for the matrix Green's function (GF) 
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E{Ko + A)f 3 ] 5 a - (r 2 + (SbiSbf)) t a f 3 - M u (u)} < ^|^+ > w = <%f . 
i 



( 



Here f 



101 . / 1 \ a (10 

T\ = \ , f 3 = I are Pauli matrices and Mn denotes the 



V 



o l J V 10 / \o -l 

matrix self-energy. Due to complicated interactions in the Hamiltonian Mu contains a number 
of terms. Here we write down few most important contributions. 
The contribution from fermion-phonon scattering is given by 

^V) = E E + V? s 5 tr )(r 2 T«' ]S , + V$6 VJ )n « t#M«?V# > w r 3 . (9) 

i'j' aa'ss' 

The contribution from linear boson-fermion scattering reads 

= E rH^h < (f>f f ip f \(f>^ f 3) (10) 



with 0g = 56; + Sbf. 

Related, but probably less important contribution from quadratic boson-fermion scattering 
(from H") 

M-f(u) = J2 Uj'U'jh < 8bi5b^ r \5b v 5b^^ > w f 3 . (11) 

This contribution is absent if one approximates H" in mean field like manner. 

There are also two contributions from three particle scattering events: fermion-phonon- 
slave boson. From H' part of the Hamiltonian we get 

Mjf"» = E E rX%Tf jsl r 3 « 0g,<^|0^-<;^ » w f a , (12) 

i'j' aa'ss' 

while from second term in H" 

i'j' aa'ss' 

To get real frequency axis Eliashberg equations in a standard formic one has to express the 
higher order Green's function appearing on the rhs' of various selfenergy pieces by the GF's 
of fermions <C ip%\ipt Xi bosons 5b^\5bt and ,,phonons" <C <j)$v\<t>-& This is 

easily achieved with help of spectral representation and decoupling of various time correlation 
functions as e.g. 

(^t)M*)M*W(o)Mo)V^ 
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This approximation neglets vertex corrections. In the present situations the vertex corrections 
stem from fermion-phonon, fermion-boson and also more complicated interactions. As already 
mentioned the assumption of small electron-phonon coupling allows us to neglect fermion- 
phonon vertex corrections. The same is not true for the fermion-boson vertices. They do 
not contain small parameters and should be taken into account. Their importance has also 
been stressed in previous studiesi'S. We shall postpone discussion of vertex corrections and 
concentrate on the anisotropy of Eliashberg equations. 

The next step is Fourier transform of the various selfenergy parts. Defining the Fourier 
components of Mij by 

s» = ^£ e ^ M) ^>)> (13) 

ij 

we get 

k'qv 1 3 

^) = i E ( e L/3 « ». r 3 i £ e^'+-W'-^) 



N 3 

k'q,p 

1 



— £ e i(fe+P-<£)(iV-ity) I _ ( 14 ) 
^ i'i' J 

Here = a gi/ + a+ gi , denotes the phonon field, a+ is creation operator of the phonon. 
Electron-phonon matrix elements and M^ qp are defined as 

a 

(15) 

mj 2 = vy\gM\ 2 \T.^mn-q-p-nj\^ 

a 

and g v {q) = {%/ /2uj v {q) ■ M) 1 ^ 2 , v% = |Jr§, e"(g) is the a-th component of the phonon 
polarisation operator, while V~ is the Fourier transform of the ionic part of electron-phonon 
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interaction, u v (q) denotes the phonon dispersion and M — the ionic mass. V c (V$) denotes the 
strength of the covalent (ionic) part of electron - phonon interaction. 
The factors of the type 

are to be treated very carefully. As already mentioned the sum goes over all sites i' and f if we 
evaluate the normal contributions to the self-energy matrix i.e. components (1,1) and (2,2) of 
the self-energy matrix which are proportional to the GF of the type <C Oci/ a \0'cp a ^> (0,0' 
are arbitrary operators). In such a case the above sum reduces to N5^, . When calculating 
off-diagonal elements of the selfenergy which depends on the GF's like <C Oc^O'c^ 3> it is 
reasonable to assume i' and f to be nearest neighbours, as in the U — oo limit only nearest 
neigbhour pairs will probably survive. In this situation the sum reduces to 'y(k) = —e(k)//2t. 

The approximation of the sum to (i',f) being nearest neighbour sites relies on the experi- 
mental fact that the superconducting coherence length is small: of order of lattice spacing. In 
an effective theory like the one presented here it translates to nearest neighbour pairs. 

To proceed we expand matrix selfenergy Sfc(u;) as 

t k (uj) = u[l - Z k (uj)]f + fc (w)fi + Xk(u)f 3 (16) 



and write down the equations for various parts of it. We get 

th^ + cth^ 1 



1 r r th^^l _|_ rth^L 1 



1 N T _u 1 Z^(u; l ) 



( ) Im 



7T 



Vk v ; 2 J J u - ui - uo N ^ kk ' K ' 



k' 



f _lw r 2 6 g/ - ^ + A + E q -e n ,_ q ASb^b+) + x«M 

where we denoted D^fa) = [^i%(^i)] 2 -[0g,(^i)] 2 -[r 2 e^-//+A+Eq ejfc/_^M & J)+Xfc'(^i)] 2 
and 
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= E \M^ p \ 2 (-hlm « ^_*r, ». 2+io , (18) 

k tk' M = Jf E 2 ("^ /m « ^>'^ 7(* - 9 - (19) 

Completely analogous set of equations is obtained for contributions to E-?(u;). The only differ- 
ence is that Kf- and Rf- will be replaced by K? P and iff given by 

= 4(-") /m « + 4(-^) /m < I >- 2 +-, (20) 

and 



= ^ E e L/--) J ™ « * W 7 (£ - q ~ k>) 



N Y 9 tt' 



E(~) ^«56t|%» W2+i0 ^k + q-k'). (21) 



Due to complicated structure of the GF entering the expression for M| A B (c<j) the corresponding 
Eliashberg equations on the real frequency axis have more complicated form 

If, f J f J th^(cth^ + cth^) + cth^cth^ + l 



uj[1 - Z^ h B {uj)] = - J dui J du 2 J du 3 - 



ivh _ Bl s l fi [■ [■ th^(cth^ + cth^) + cth^cth^ + 1 
0f "(u) = - dwi duj 2 / dco 3 2 — 2 2 — 2 - 2 

k A J J J UJ — UJ] — UJo — UJ?, 



vh _ Bl , If- [■ [■ th^(cth^ + cthSf.) + cth^cth^ + 1 

Xf B (UJ) = - / dw l / / dw 3 — — 

fc A J J J UJ — UJ\ — uj 2 — ^3 



^E^ K^)(--)/m ^ , (22) 



with 



Kg,-> 2 , a*) = w £ {r 2 |M^_/(--) Jm « fe -j0_^ » W2+M (-I)/m « %|<ft± ». 3+ *o 
+ r2 |^, w,ol 2 (-^) J ™ « ^J*-iE-*„ ("^) Im « ^1% »-3 + ,o} (23) 



and 



Kff B (uj 2 , ^ 3 ) = ^ E /m « (-^) /m « * W »«-»+* • 

•7(fc - g) +r 2 |M|_ =0 | 2 (-i)/m < >^ 2+i0 < 56t|%> W3+io 7 (£ + p- g)| 



(24) 
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Equations (17-24) supplemented with expressions for phonon and slave - boson Greens func- 
tions form a complete set and have to be solved selfconsistently. 

There are few differences between them and usual Eliashberg equations for superconductors. 
The most important is the presence of two different types of kernels Kkk' and Kkk 1 determining 
normal {Zy,Xk') and anomalous (0^0 parts of the self-energy, respectively. 

Other point worth to mention is strong renormalisation of the bare electron spectrum eg = 
—2t^(k). At the mean field level of the theory it is both shifted and mass renormalised to 
r 2 e£ + A. Beyond mean field level for bosons the nontrivial energy dependence comes from 
renormalisation of the fermion (matrix) Greens function given by the Xg ( a; ) factor. Note, that 
from the equation for xf (^) an d definition (20) of K?^, (u) we have 

X fH = CH+e|CxH, (25) 

where C(u) and C\(uj) are weakly k dependent functions to be calculated from equation (17). 
To take this effect into account one has to solve the equations at each point of the Brillouin 
zone — a very difficult task. 



IV. RESULTS AND DISCUSSION 



The complete solution of Eliashberg equations which are integral equations with complicated 
kernels would require finding all functions Z k (u), 0fc(cu), XfcO^O at all frequencies and all points 
in the Brillouin zone. In the standard approach the valid argument was that most important 
changes are to be expected at the Fermi level. The various pieces of selfenergies were thus 
averaged over the wave vectors lying on the Fermi surface. In this way it was possible to reduce 
the equations to the integral equation of the single variable: the frequency. It is not obvious 
that the same arguments do apply in the present case especially to the boson contribution to 
self-energy. 

Sticking to the mean field approximation for slave bosons leads to 5b q = 5b+ = 0. In this 
approximation only phonon contribution to Eg survives, Sfc(c^) = Ti^ h (uj). In this case most of 
arguments developed previously^ apply and one can use known parametrisations of solutions 
to Eliashberg equations. The relevant parameters describing superconducting transition tem- 
perature T c are the typical phonon frequency ujo and electron-phonon parameter A e _ p h- The 
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last parameter is determined by the kernel K ph . In our theory we do have two different kernels 
as also found previously by Zieli/nski and coworkers^. 

We thus define two coupling constants A c - P h and A e _ p h- The first one is related to wave 
function renormalisation and other describes genuine superconducting coupling in strongly 
correlated system. We define A e i_ p h by0 

A c -ph = 2 I — a 2 F{uj) (26) 
J UJ 

where 

°? F (") = E K H H S(r 2 e n -fi + A)5(r 2 e f : l - p + A) // £ S(r\ - p + A) 
kP k 

Similar expression defines A e - P h except that K% k ,(u) enters. In the general case there exist four 
more parameters: A e _ B , \ £ -b, K-ph-B and \ e _ ph _ B . 

For numerical purposes it is convenient to take optical phonons with uj qv — u . In this case 
the frequency integral in (26) can be calculated very easily. 

Strong wave vector dependence of the kernels Kg^, makes the question of symmetry of the 
order parameter very important. To get some information on this point we calculate various 
symmetry components A*_ p/l of the pairing interaction. We define them as 

K- Ph = ^ E Kp(?\ - /* + A) 5{r% - f, + X)g i (k)g i (k')//A i 

kk' 

We have taken gi = cos k x a + cos k y a, g<i = cos k x a — cos k y a, g% = sin k x a + sin k y a and 
g 4 = sin k x a - smk y a; At = T,k9i{k)S{r 2 e l : - fi + A). 

There are few parameters which control the behaviour of the system. These are: carrier 
concentration n, strength of the electron - phonon interaction parameters V c and V{, the fre- 
quency ujq of the optical phonon etc. All the energies are measured in units of D = 2t, and for 
numerical purposes we have taken ujq = 0.02, and the value of the ionic mass M appropriate 
for oxygen. 

In this work we shall concentrate on the results obtained from the kernel , describing 

kkr ° 

anomalous part of the self-energy. In the previous work! we have analysed in some detail the 
properties of the system resulting from the kernel entering normal self-energy. 

We start the discussion by showing the results obtained for Vi = 0. The covalent part 
of the interaction is strongly k dependent. It is described by the matrix element |M^| 2 in 
equation (15). In figure (la) we show the results obtained by making the same approximation 
as in previous treatment of the similar model,@ where |M^| 2 has been replaced by wave vector 
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independent constant multiplied by r 4 . The concentration dependence of s, p and d components 
of the coupling constant are shown by solid, dotted and dashed lines, respectively. It is the 
factor r 4 = (1 — n) 2 , which strongly supresses the values of A; at higher electron concentrations. 
The important point is that for small doping of holes the d-wave component is the largest one. 
This signals the possibility of the d-wave pairing in this parameter range. For 5 = 1 — n > 0.3 
the p-wave pairing takes over. At quite small electron concentrations extended s-wave seem to 
be most stable pairing. 

In figure (lb) we show similar set of data obtained with |MgJ 2 replaced by its Fermi surface 
average M\. Due to oscilatory dependence of |M^-| on wave vectors the average value of it 
is strongly concentration dependent (solid line in figure (lb)) and reduced in comparison to 1 
(the value arbitrarily assumed for data shown in Fig. (la)). Small values of M F at large n are 
due again to r 4 factor, while at small n to phase space restrictions. 

Taking the full wave vector dependence of |M^J leads to the results shown in figure (lc). 
Solid line in this figure shows total coupling A c i_ p h calculated from (26). Other lines represent 
the symmetry components \ as indicated. The values of all couplings are very small. Thus 
in the mean field approximation for slave bosons one is not able to get realistic values of the 
superconducting transition temperature from covalent part of the electron-phonon interaction. 
Let us look at the relative values of the various couplings. At large electron concentration 
(small hole dopping) the d- wave component of A exceeds by far all other components. 

Bigger values of coupling constants are obtained when ionic part of the electron - phonon 
interaction is taken into account. Figure (2a) shows the concentration dependence of A (solid 
line) and \ (dashed lines) obtained for V c = and Vi = 1 (in units of D). The values of the 
coupling parameters are quite realistic now. The nice feature of this result is the possibility of 
d symmetry pairing at small hole concentration. The presence of both ionic and covalent parts 
of electron-phonon interaction leads to small changes of the phase diagram (cf. Fig.2(b)) . The 
d symmetry remains most probable at small hole dopings 5 < 0.3. For intermediate carrier 
concentrations p-wave pairing is most stable, while extended s-wave is stable at small electron 
concentration n < 0.3. The slight decrease of the values of couplings as compared to figure (2a) 
is due to interference effects between interaction channels (cf. expression (15) for |M^| 2 ). In 
this work we have concentrated on the calculation of the electron - phonon coupling constants. 
To obtain realistic estimations of the concentration dependence of the superconducting transi- 
tion temperature one has to take simultaneously into account the condensation of the (slave) 
bosons. Boson condensation temperature Tf of the 3d massive bosons is proportional to the 
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2/3 power of the boson concentration (1 — n). For the masless bosons the power changes to 
1/3 in 3d case and 1/2 for interesting case of two dimensional system. Becouse superconduct- 
ing transition temperature, in the systems with covalent coupling (c.f. figures la,lb,lc), also 
vanishes for n — > 1, the additional condition for Bose - Einstein condensation can change the 
phase diagram in quantitative way only. The situation is different in systems with ionic or 
mixed (covalent plus ionic) interactions and especially for the phase with ci-wave symmetry of 
the order parameter. The coupling constant takes on sizeable values (c.f. figs. 2a and 2b) 
at low hole concentration. In this case the additional condition will modify the T c vs.n phase 
diagram making it qualitatively more similar to that obtained for covalent electron - phonon 
interaction. 

In conclusion we have presented the derivation of the Eliashberg equations for the system 
with strong electron - electron and electron - phonon interactions. The results show the impor- 
tance of the correlations which make the superconducting order parameter very anisotropic. Its 
symmetry depends on the carrier concentration. We have illustrated the theory by showing the 
results of calculations of coupling constants in the mean field approximation for slave bosons. 

As already mentioned beyond mean field level there are important changes of the structure of 
the theory. The mixed fermion - phonon - slave boson contribution to the selfenergy and corre- 
sponding coupling constants will overtake the fermion - phonon contributions considered here. 
The x component of the selfenergy due to fermion - slave boson scattering acquires nontrivial 
energy eg dependence. The study of both these effects will be reported in the forthcomming 
paper. 
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FIGURES 

FIG. 1. The concentration dependence of the coupling constants Aj calculated for a model with 
Vi = 0, V c = 1 and: (a) constant value of the electron - phonon matrix element \M^-\ = 1.0, (b) \M^-\ 
averaged over the Fermi surface and (c) the actual value \M^~\ with full wave vector dependence taken 
into account 

FIG. 2. The pairing coupling constants vs n for a model with: (a) ionic electron - phonon inter- 
action (Vi = 1, V c = 0) and (b) both ionic and covalent interactions (Vi = 1, V c = 1) 



17 



